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Geophysical applications of the Global Positioning System (GPS) require the ca- 
pability to estimate and propagate satellite orbits with high precision. An accurate 
model of all the forces acting on a satellite is an essentia/ part of achieving high 
orbit accuracy. Methods of analyzing the perturbation due to thermal radiation 
and determining its effects on the long-term orbital behavior of GPS sate//ites are 
presented. The thermal imbalance force , a nongravitational orbit perturbation pre- 
viously considered negligible, is the focus of this article. The Earths shadowing of 
a satellite in orbit causes periodic changes in the satellite’s thermal environment. 
Simulations show that neglecting thermal imbalance in the satellite force model 
gives orbit errors larger than 10 m over several days for eclipsing satellites. This 
orbit mismodeling can limit accuracy in orbit determination and in estimation of 
baselines used for geophysical applications. 


I. Introduction 

Mismodeling of satellite force parameters can have a 
significant effect on satellite orbits, especially in orbit pre- 
diction [1], Some applications require the capability to 
estimate and propagate satellite orbits with high preci- 
sion. TOPEX/POSEIDON precision orbit determination, 
for example, requires precise modeling of nongravitational 
forces to fulfill mission requirements [2]. In addition, some 
of the observed drag and orbit decay on the Laser Geo- 
dynamics Satellite (LAGEOS) spacecraft have been at- 
tributed to unmodeled thermal forces [3,4]. To achieve 
a high level of orbit accuracy, an accurate model of all the 
forces acting on an Earth-orbiting satellite is necessary. 

The focus of this analysis was to assess the effects 
on satellite orbits of neglecting thermal reradiation and 


of mismodeling nongravitational forces. Radiative heat 
transfer between a satellite and its environment is the ba- 
sis for the thermal force model. A satellite in Earth orbit 
is continuously illuminated by radiation, most of which 
comes from the sun. The thermal imbalance force is di- 
rectly related to the temperature distribution of the satel- 
lite in its changing environment. An uneven temperature 
distribution causes surfaces to reradiate energy at different 
rates. Some studies have shown that most of the thermal 
gradient forces on the TOPEX/POSEIDON satellite orig- 
inate within the spacecraft body [2]. Other analyses have 
shown that the dominant source for thermal reradiation 
forces on a Global Positioning System (GPS)-like satellite 
is the solar panels, due to their large exposed area and 
low heat capacity [5]. The satellite’s heated body rera- 
diates energy at a rate proportional to its temperature, 
losing the energy in the form of photons. By conservation 
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of momentum, a net momentum flux out of the body cre- 
ates a reaction force against the radiating surface, and the 
net thermal force can be observed as a small perturbation 
that affects long-term orbital behavior of the spacecraft 
[5]. The partial differential equations and boundary con- 
ditions describing the temperature distribution and the 
heat transfer between surfaces, along with the application 
of the finite element method, are presented in this article. 
A brief description of the statistical estimation technique 
used for studying the effect of the thermal imbalance force 
on satellite orbits is included. 


II. Radiation and Heat Conduction 
Formulation 

Two types of heat transfer that affect a spacecraft’s 
orbit are radiation and heat conduction. The exchange 
of energy between the spacecraft and its surroundings 
is described by radiation heat transfer. Conduction is 
the transfer of heat by molecular motion within a solid 
medium. Figure 1 shows these types of heat transfer. 
The rate of radiant energy transfer is given by the Stefan- 
Boltzmann Law [6]: 

E r = eaT* (1) 

By conservation of momentum, the thermal force or rate 
of change of momentum for a radiating surface element, 
assuming a Lambertian surface, is expressed as [5] 

2 eaT* f v 

dfthermal — 77 dAtl a (2) 

o C 

The unit vector in Eq. (2) is defined as normal to and 
directed out of the sun-tracking surface of the solar panel. 
The differential force must be integrated over the entire 
surface to determine the complete thermal force: 

fthennal = J eT? dAli a (3) 

Clearly, thermal forces cannot be computed unless space- 
craft surface temperatures are known. In general, the tem- 
perature at any point within a body satisfies the heat equa- 
tion [5] 

KV 2 T= P C p ^ (4) 


The solution to this second-order partial differential 
equation requires that boundary conditions be specified. 
The boundary conditions are defined by thermal radiation 
and heat conduction. As given by the conservation of en- 
ergy principle, the total amount of energy coming into a 
surface is equal to the total amount of energy leaving the 
surface, assuming there is no internally generated or stored 
energy (no sinks and no sources). The boundary condition 
for the satellite surface can be obtained by expressing this 
condition as 

9in = <7out (®) 

where q in is the amount of incoming radiative energy due 
to external sources and internal conduction, and q out is the 
amount of radiative energy leaving the boundary due to 
reradiation and conduction. Figure 2 shows the conserva- 
tion of energy principle for a satellite solar panel surface. 
Using this concept, the boundary conditions for each sur- 
face were constructed. The incident radiative solar energy 
received per unit area per unit time by side a and side b 
of the solar panel are represented by h a and h\, [5]. 


KAj± = etvAT? - h b A 

(6a) 

- KA = e a <rATt - h a A 

(6b) 


The actual amount of incident radiative energy received 
by each side of the solar panel is a function of panel orien- 
tation and the orbit of the satellite. The subscript a rep- 
resents the left boundary in local coordinates (cold side), 
and the subscript b represents the right boundary, which, 
for a GPS-type satellite, is assumed to be continuously 
facing the sun during orbit. The term on the left side of 
the equal sign in Eqs. (6a) and (6b) is the heat flux, en- 
ergy per unit time per unit area, in the local x-direction, 
which is perpendicular to the solar panel face. The values 
used for some of the parameters described above are listed 
in Table 1 and are consistent with values used for GPS 
satellites. 

PDE-Protran, a finite element method program, was 
used to solve the transient heat conduction and radiation 
problem presented here. PDE-Protran was developed by 
Granville Sewell and is a general-purpose two-dimensional 
partial differential equation solver [7]. This software was 
combined with a program which incorporated material 
properties, the satellite’s orbit orientation, and thermal en- 
vironment to determine solar panel surface temperatures. 
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Grid points were chosen to divide the solar panel into small 
sections or “elements,” where the temperature of the so- 
lar panel was computed for each grid point in one dimen- 
sion, across the thickness of the solar panel. These grid 
points coincide with the boundaries between each layer 
of the solar panel’s “sandwiched” materials (listed in Ta- 
ble 2). Accurate and current knowledge of physical pa- 
rameters such as surface emissivity, thermal conductivity, 
heat capacity, and material density is required. For this 
analysis, the material properties are assumed to remain 
constant throughout the satellite’s orbit, and only the so- 
lar radiation environment varies with time as the satel- 
lite experiences eclipsing or shadowing from the sun by 
the Earth. The material parameters directly influence the 
thermal forces which are calculated and have an effect on 
the prediction and propagation of the spacecraft trajec- 
tory. Also, these material properties may change in time 
or degrade, due to the harsh environment of space. 1 

III. Orbit Analysis Technique 

In this investigation, the equations of motion for an 
Earth satellite are assumed to include the two-body grav- 
itational effect and the thermal imbalance forces only, and 
are given in vector form by 

£*r ^thermal /«\ 

r = =■ H IG 

r d m 

and thermal imbalance force perturbing the satellite, 
ft hermai j i® computed as 

- 2 a A 

^thermal = — (ol? ~ C«7?) U a (8) 

The effect of the thermal imbalance force on a satellite 
can be observed by comparing the perturbed orbit with 
the unperturbed two-body orbit in time. Since there is no 
closed-form analytical solution for the perturbed equations 
of motion, a numerical integration technique was neces- 
sary to solve the ordinary differential equations of motion, 
Because the perturbed and unperturbed orbits originate 
with the same initial conditions, the displacement between 
them at a given time can be observed. A least-squares es- 
timation technique is used to determine the state of the 
satellite in its orbit at a specified epoch [8]. The initial 
conditions of one orbit can be adjusted at a given time 
to eliminate the secular divergence between the perturbed 


1 T. T. Lam, personal communication, Aerospace Corporation, El 

Segundo, California, May 1990. 


and unperturbed orbits to observe the periodic behavior. 
In this analysis, two types of GPS satellite orbits were 
studied. The satellites of the GPS are distributed in six, 
evenly spaced orbit planes. When completed, the final con- 
stellation will consist of 24 satellites at an orbit altitude of 
approximately 20,000 km, with an orbit period of about 
12 hr. In this constellation, most satellites are exposed 
to full sunlight. As the orbit geometry changes, however, 
some GPS satellites will experience eclipsing or shadowing 
from the sun by the Earth. Both eclipsing and noneclips- 
ing satellites are the focus of this study. Throughout its 
orbit, the GPS solar panel maintains a fixed orientation 
toward the sun. Nodal motion was not considered, since it 
is not significant for the short time interval (1 week) used 
in this study. No internally generated energy was mod- 
eled in this study, but the absorbed solar radiation that is 
converted to electricity was modeled, using the efficiency 
of the solar panel at 14.1 percent. Although studies have 
shown that for a TOPEX/POSEIDON satellite the ther- 
mal radiation forces originating with the spacecraft body 
are twice those from the spacecraft solar panels, the major 
source for thermal reradiation forces on a GPS-like satel- 
lite is the spacecraft’s thin, large, solar panels [2,5]. Con- 
sequently, in this analysis, the GPS satellite’s main body 
was not considered. Other studies are currently consid- 
ering this problem of modeling thermal reradiation forces 
for a complete GPS spacecraft. 

IV. Discussion of Results 

In order to determine the direction and magnitude of 
the thermal force, the surface temperatures were calcu- 
lated using the finite element method program, PDE- 
Protran [7]. Several simulations were tested. The input 
required for the simulation is shown in Table 2. This table 
lists the material properties for a Block II GPS satellite so- 
lar panel [9,11]. 2 “ 5 The initial conditions included a solar 
panel orientation perpendicular to the sun and an initial 
temperature of 300 K. The time step used in the analysis 
was 100 sec (one GPS orbit is approximately 43,200 sec, 
and the eclipsing period lasts approximately 3200 sec). 

GPS satellites experience an eclipsing season for only 
a few weeks every year. Eclipsing has a strong effect on 
the solar radiation environment of those satellites. This is 


2 Ibid. 

3 W. Pence, personal communication, Rockwell International, Seal 
Beach, C alifomia, October 1990. 

4 J. Albeck, personal communication, Spectrolab Corporation, Syl- 
mar, California, August 1990. 

s D. Marvin, personal communication, Aerospace Corporation, El 
Segundo, California, January 1991. 
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evident in the temperature of a GPS satellite solar panel 
over one orbit, as shown in Fig. 3. The steady-state tem- 
perature for the sun- facing side is approximately 317 K, 
and for the shaded side it is 313 K. These values compare 
well with the approximate value of 313 K, which has been 
measured on the cold, shaded side of the solar panel for a 
GPS satellite. 6,7 The face exposed to the sun has not been 
directly measured and, therefore, the temperature differ- 
ence between the surfaces is not well known, but is believed 
to be approximately 5 K [10]. 8_1 ° During the eclipse pe- 
riod, which lasts approximately 1 hr, panel temperature 
declines approximately 253 K. After exiting the shadow 
region, the solar panels slowly return to their steady-state 
temperatures over approximately 3 hr. 

Modeling the cover-glass surface accurately has been 
difficult during this study since data on the properties 
of this surface were not readily available. The thermal 
conductivity of this fused silica layer is very low as com- 
pared to that of two other dominant layers, the aluminum 
core and the solar cell layers. 11 The cover-glass layer, on 
the sun-facing side of the solar panel, contributes most of 
the temperature imbalance, primarily because of its low 
thermal conductivity and greater thickness as compared 
to other solar panel layers, especially the aluminum core. 
Although it is believed that the solar panel cover-glass 
layer is transparent to all incident radiation, the material 
properties of this specific layer of the solar panel were not 
deleted from this analysis. It was important to simulate 
the solar panel as it exists in orbit to observe the long- 
term orbital effects of the thermal imbalance force on a 
GPS satellite. This simulation is adequate as long as the 
correct material properties are used in the analysis. 

As an example, two simulations were performed using 
identical solar panel parameters (values given in Table 2) 
except for different thermal conductivities for the sun- 
facing cover- glass layer. These simulations are presented 
to show the sensitivity of the temperature calculations to 
the thermal conductivity of the cover glass. The value for 
the thermal conductivity given in test case 2, shown in 
Table 3, was used to demonstrate how unrealistic thermal 
forces can be computed when using incorrect values for 
the solar panel material properties. Previously, however, 
this was believed to be the correct value for the thermal 
conductivity of the fused silica cover-glass layer of a GPS 


6 Pence, op. cit. 

7 Albeck, op. cit. 

8 Pence, op. cit. 

9 Albeck, op. cit. 

10 Marvin, op. cit. 

11 Ibid. 


satellite solar panel [9]. 12 The results, shown in Table 3, 
describe the steady-state temperatures and thermal ac- 
celerations that were computed using the specified values 
for the cover-glass thermal conductivity. Again, both test 
cases shown in Table 3 are identical except for the value of 
thermal conductivity for the solar panel cover- glass layer. 

In this article, the reference frame is defined as space- 
craft-centered radial and along-track components. The 
along-track component is also referred to as the trans- 
verse or down-track direction, defined in the direction of 
the satellite velocity vector. Figure 4 shows radial and 
along-track components of the acceleration due to thermal 
reradiation over one orbit for an eclipsing satellite. These 
components compare well with those of studies which have 
shown unmodeled nongravitational forces to cause errors 
of the magnitudes shown in Fig. 4 [10]. Also, these results 
were computed using the information presented in Table 2 
and described as test case 1 in Table 3. 

Figure 5 shows the differences between two orbits, one 
computed using two-body effects only and the other — an 
eclipsing orbit — computed with two-body effects and the 
thermal imbalance force during 1 week. The radial rms is 
0.5 m, and the along-track rms is 5.2 m. These results were 
computed using a technique similar to the method used to 
predict satellite orbits based on a set of initial conditions 
and a complete force model of the spacecraft, which could 
include the solar radiation pressure and thermal imbalance 
force. In the catse of the eclipsing satellite after 7 days, the 
along-track components differ by approximately 13 m. 

Figure 6 also shows the differences between two orbits, 
one computed using two-body effects only and the other 
computed using two-body effects and the thermal imbal- 
ance force for a satellite not in an eclipsing plane. The 
radial rms is 0.5 m, and the along-track rms is 1.6 m. It 
can be seen from these results that an eclipsing satellite 
experiences a larger perturbation in the along-track direc- 
tion over the span of 1 week than a satellite not in an 
eclipsing orbit plane. For the noneclipsing satellite after 7 
days, the along-track difference is approximately 5 m. 

Figures 7 and 8 represent the results computed using 
a least-squares estimation algorithm in which the simu- 
lated observation data contained only the two-body grav- 
itational and thermal imbalance reradiation forces. The 
force model used in the estimation algorithm contained 
the two-body gravitational force model and a solar radi- 
ation pressure model to observe the ability of the force 
model to account for thermal imbalance forces which have 


12 Lam, op. cit. 
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been difficult to model but exist in the observations. The 
calculated best estimate of the satellite epoch state, in the 
least squares sense, includes the satellite position, velocity, 
and a solar-radiation pressure scale factor. 

Figure 7 shows the orbit fit residuals for a satellite in 
an orbit plane that is regularly eclipsing. The radial rms 
is 5 cm, and the along-track rms is 80 cm. After 7 days, 
the along-track orbit error is almost 2 m. These results 
show that the solar- radiation pressure scale factor in the 
estimation scheme is capable of absorbing most of the orbit 
error due to thermal reradiation, but not all of the orbit 
error, especially in the along-track direction. 

Figure 8 also shows orbit fit residuals for a GPS satellite 
using the same estimation technique, but for a satellite in 
a noneclipsing orbit plane. The radial rms is 9 mm, and 
the along-track rms is 17 cm. After 7 days, the along-track 
orbit error is approximately 40 cm. Clearly, the eclipsing 
of the satellites has an influence on the orbit errors when a 
thermal reradiation force is not included in the estimation 
force model. Larger orbit errors are calculated when the 
satellite is in an eclipsing orbit plane. A 1-week prediction 
can be made using the satellite state computed for the best 
least squares estimate in Fig. 7 and compared to the best 
least squares estimate for that predicted week. Studies 
have shown that, for eclipsing satellites, the quadratic-like 
growth in the along-track direction can give errors as large 
as 50 m after a 1-week prediction [12]. 

V. Concluding Remarks 

The current analysis has shown that orbit errors larger 
than 10 m occur when mismodeling nongravitational forces 
such as the thermal imbalance force presented here. A 
finite element method technique has been used to calcu- 
late satellite solar panel temperatures which are used to 
determine the magnitude and direction of the thermal im- 
balance force. Although this force may not be responsi- 
ble for all of the force mismodeling, conditions may work 
in combination with the thermal imbalance force to pro- 
duce accelerations on the order of 1.0 x 10 9 m/sec 2 . One 
possible contribution currently being studied is the solar 


panel misalignment acting together with the thermal im- 
balance force, a contribution that may account for much of 
the unmodeled perturbations. If submeter-accurate orbits 
and centimeter-level accuracy for geophysical applications 
are desired, a time-dependent model of the thermal imbal- 
ance force should be used, especially when satellites are 
eclipsing and the observed errors are larger than those for 
satellites in noneclipsing orbits. One study has shown that 
estimating additional stochastic solar radiation parame- 
ters improves GPS orbit accuracy significantly, especially 
for eclipsing satellites [13]. This technique can be used 
to absorb the orbit error caused by mismodeling thermal 
imbalance forces. 

Although modeling the spacecraft solar panels alone 
may be considered insufficient, thermal force modeling of 
the entire spacecraft is a complicated problem. This has 
been done for spacecraft such as TOPEX/POSEIDON, 
where precise orbit determination is critical to mission suc- 
cess [2]. The study presented here, however, focused only 
on modeling the solar panels, where the material composi- 
tion is not nearly as complex. Also, the problem of radia- 
tion absorbed and conducted through the solar panel and 
reradiated out is a simple one-dimensional time-dependent 
heat- transfer problem, with no internal heat generation 
from scientific instruments or electronics. 

Nongravitational perturbations like the thermal imbal- 
ance force have been observed for years on satellites like 
LAGEOS, and are still not completely understood. Ther- 
mal forces are dependent on the environment and specif- 
ically on such parameters as the satellite mass, cross- 
sectional area, and material composition. Unfortunately, 
these parameters can change or degrade with long-term 
exposure in space. For this reason, it may be more appro- 
priate to estimate stochastic force parameters to represent 
the thermal reradiation forces since the nature and rate of 
material degradation of the satellite in orbit is unknown 
[13]. The results obtained using the finite element model 
in this study agree with the work of others who have con- 
ducted similar studies using the finite difference technique 
to determine spacecraft thermal gradient forces in an effort 
to improve the satellite force models. 
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Appendix 

Glossary 


Character 

A 

c 

C p 

E r 

f 

h 

K 

m 

n 

<1 

T 

r 

t 

T 

a 

£ 

P 

* 

<T 

P 

Subscript 

a 

b 

in 

out 

r 

thermal 


Purpose and/or usage 
surface area 
speed of light 
specific heat 

energy emitted by a real body, summed over all wavelengths 

thermal imbalance force per unit area 

incident solar radiation received by solar panel 

thermal conductivity 

satellite mass, kg 

unit vector normal to surface of solar panel 
radiative energy 

geocentric satellite position vector 

radial distance from Earth’s center of mass to satellite 

time 

temperature 

surface absorptivity 

emissivity 

density 

solar constant 

Stefan-Boltzmann constant 

Earth’s gravitational parameter 

solar panel sun-tracking front side 
solar panel back side, no direct sunlight 
incoming to the spacecraft 
leaving the body 
radiative 

thermal imbalance 
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Table 1. GPS thermal and orbit parameters. 


Model parameter 

Value 

Initial orbit radius, m 

26,550,000 

Surface emissivity (e a ) 

0.78 

Surface emissivity (ej,) 

0.83 

Surface absorptivity (<*a), percent 

0.77-14.1 (panel efficiency) 

Solar panel surface area (A), m 2 

10.832 

Satellite mass (to), kg 

845 

Ini tied panel temperature (t = 0), K 

300 

Stefan -Boltzmann constant (<t), W/m 2 K 

5.6699 x 10” 8 

Speed of light (c), m /sec 

2.998 X 10 8 

Solar constant ('I'), W/m 2 

1368.2 

Total panel thickness (eight layers), m (in.) 

0.01478 (0.582) 


Table 2. GPS Block II solar panel properties. 


Panel layer 
composition 

Thickness, 

m 

Density, 

kg/m 3 

Specific heat, 
J/kgK 

Conductivity, 

W/mK 

Cover glass 

0.00749 

2186.622 

753.624 

1.417 

Adhesive 

0.00005 

1079.472 

1256.04 

0.116 

Solar cell 

0.00025 

2684.84 

711.756 

147.994 

Interconnect cell adhesive 

0.00018 

1051.793 

1256.04 

0.116 

Kapton cocured 

0.000076 

1162.508 

1130.436 

0.1506 

Graphite epoxy 

0.00019 

2186.622 

1373.27 

0.8706 

Aluminum core 

0.00635 

24.91088 

1046.7 

250.966 

Graphite epoxy 

0.00019 

2186.622 

1373.27 

0.8706 


Tabla 3. Cover glass thermal simulation. 

Test case 1 

Test case 2 

K = 1.417 W/mK 
Hot side T a = 317.41 K 
Cold side T \ = 313.66 K 

Thermal acceleration = 1.88 X 10“ 10 m/sec 2 

K = 0.04327 W/mK 
Hot side T a = 340.30 K 
Cold side T h = 285.37 K 
Thermal acceleration ~ —8.01 x 10 -9 m/sec 2 
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Fig. 3. Temperature history simulation for a GPS solar panel. 


Fig. 1. General heat transfer diagram for a spacecraft 



Fig. 2. Conservation of energy diagram for a satellite surface. 



Fig. 4. Radial and along-track components for the thermal force 
over one orbit 



Fig. 5. Radial and along-track orbit differences, eclipsing satellite. 
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Fig. 6. Radial and along-track orbit difference#, noneclipsing satellite. 



TIME, days 

Fig. 7. Orbit fit residual*, with a solar-radiation pressure scale 
factor, eclipsing satellite. 



TIME, days 


Rg. 8. Orbit fit residuals, with a solar-radiation pressure scale 
factor, noneclipsing satellite. 
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